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ABSTRACT 

Let u(x,t) be a weak solution of the Euler equations, governing the 

inviscid polytropic gas dynamics; in addition, u(x,t) is assumed to respect 

the usual entropy conditions connected with the conservative Euler 

equations. We show that such entropy solutions of the gas dynamics equations 

satisfy a minimum entropy principle , namely, that the spatial minimum of their 

specific entropy. Ess inf s(u(!c,t)), is an increasing function of time. This 

x 

principle equally applies to discrete approximations of the Euler equations 
such as the Godunov-type and Lax-Fried richs schemes. Our derivation of this 
minimum principle makes use of the fact that there is a family of generalized 
entropy functions connected with the conservative Euler equations. 
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1 . INTRODUCTION 


Many phenomena in continuum mechanics are modeled by hyperbolic systems of 
conservation laws 

a, d 

gPjT + l — = o, (x = (x 1 ,x d ),t)e Rx[0,«), (1.1) 

k=l x k 

where are smooth nonlinear flux mappings 

T 

of the N-vector of conservative variables u= u(x,t) *= (u^,**«,u^) . 
Friedrichs and Lax [3] have observed that the hyperbolic nature of such models 
is revealed by the property of most of those systems being endowed with a 
generalized 


Entropy Function ; A smooth convex mapping U(u) augumented with entropy flux 
mappings F = F(u) = (F^ ^ (u) , • • * ,F^^ (u) ) , such that the following 
compatibility relations hold 


yT o p(k)T 

u u u 9 


k = 1 ,2 , • • • ,d. 


( 1 . 2 ) 


T 

Multiplying (1.1) by and employing (1.2), one arrives at an equivalent 
formulation of the compatibility relations (1.2), namely, that under the 
smooth regime we have on top of (1.1) the additional conservation of entropy 




ar< k > 

3x k 


0 . 


(1.3) 


(k) 

Owing to the nonlinearity of the fluxes f (u) , solutions of (1.1) may 
develop singularities at a finite time after which one must admit weak 
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solutions, i.e., those 
conservative equations, 
regularized problem, 


derived directly from the underlying integral 
Considering (1.1) as a strong limit of the 


. d 

|^+ y 
9t w 


3 f 


(k) 


? 3 2 u 

V l — T » 
k=l 3x, 


y + 0, 


(1.4) 

y 


then following Lax [9] and Krushkov [8], we postulate as an admissibility 
criterion for such limit solutions an entropy stability condition which 
manifests itself in terms of an 


Entropy Inequality ; We have, in the sense of distributions, 


8U + j 3F (k) < 
9t k=l 9x k 


(1.5) 


Weak solutions of (1*1), which in addition satisfy the inequality (1.5) 
f° r all entropy pairs (U,? 1 ) connected with that system, are called entropy 
solutions . Having a (weakly) nonpositive quantity on the L.H.S. of (1.5) 
is thus a consequence of viewing these entropy solutions as limits of 
vanishing dissipativity mechanisms. In particular, the inequality (1.5) 
implies that the total entropy in the domain decreases in time (we assume 
entropy outflux through the boundaries) 

/_ U(u(x,t))dx _< 0. (1.6) 

x 

( ^Krushkov [8, p. 241] has termed such solutions simply as generalized 
solutions. 
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In this paper, we consider entropy solutions, 

u = (p ,*,E)^ (1.7a) 

of the Euler equations. These equations govern the inviscid polytropic gas 
dynamics , asserting the conservation of the density p , the momentum 

T m 

m ■ (m^n^m^) » and energy E. Let q = — denote the velocity field of 

such motion. Then, expressed in terms of the pressure, p, 

2 

p = (Y“l)» [e - I /2 *p |q| ], y = adiabatic exponent, (1.7b) 

the corresponding fluxes in this case are given by^' 

= (m k ,q k ‘B + p»e^ k \q k (E + p)) T , k = 1,2,3. (1.7c) 

The main result of this paper asserts that entropy solutions of Euler 
equations satisfy the following 

Minimum Principle: Let u = u(x,t) be an entropy solution of the gas 

dynamics equations (1.7) and let 

S(x,t) s S(u(x, t) ) - ln(pp”^) (1.8) 


( 2 ) 


With e 


(k) 


denoting the unit Cartesian vectors e 


(k) 


kj * 
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denote the specific entropy of such solution. Then the following estimate 
holds 

Ess inf S(x,t) > Ess inf S(x,t = 0). (1.9) 

|x|<R I x I <R+t • q 

Here q max stands for the maximal speed |q| in the domain . 

The proof of this assertion is provided in Section 3 below. Prior to that 

we elaborate in Section 2 on the entropy inequality connected with the gas 

dynamics equations. In particular, Harten [5] has shown that there exists a 

whole family of entropy pairs associated with these equations, a fact which is 

essential in our derivation of the minimum principle. 

As an immediate consequence of the minimum principle, we conclude that 

Ess inf_S(x,t) is an increasing function of t for every entropy solution of 
x 

(1.7). The following argument sheds additional light on this conclusion in 

the case of a piecewise-smooth flow. To this end, an arbitrary particle 

currently located at (x,t) is traced backwards in time into its initial 

position at t « 0. Since the specific entropy of such particle remains 

constant along the particle path — except for its decrease when crossing 

backwards shock waves, it follows that its value S(x, t) is greater or equal 

than that of the initial spatial minimum Ess inf_S(x,t = 0), as asserted. 

x 

In contrast to the above ^Lagrangian' argument, the derivation of the minimum 
principle outlined below, is purely an 'Eulerian' one. It enables us to relax 
the regularity assumption on the flow, and—since we do not follow the 
characteristics, it equally applies to discrete approximations of the Euler 


equations 
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In Section 4 we consider approximate solutions of the Euler 

equations, w(x v> t), which respect the entropy decrease estimate (1.6), 

T U(w(x ,t + At))Ax < T U(w(x ,t))Ax . (1.10) 

L y) x v v ’ V — V V 

We note that such approximate solutions are obtained by entropy stable 
schemes satisfying the cell entropy inequality 

U(v(7 ,t + At)) < U(w(x v ,t)) + l " F v-l/o ] » ( 1 * 11 > 

k“l Ax 1 /z 

v 

e.g., the Godunov-type and Lax-Friedrichs schemes [6], We have 

Minimum Principle : Let w( x^ , t ) be an approximate solution of the gas 

dynamics equations (1.8) and let 


S(x v ,t) = S(w(x v ,t)) «= lntpp^) (1.12) 

denote the specific entropy of such solution. Assume that its total entropy 
decreases in time, (1.10). Then the following estimate holds 


S(x,t + At) > Min[ S(x ,t)]. (1.13) 

V 

In the case of entropy stable schemes, (1.11), a more precise estimate is 
obtained which takes into account the support of the schemes^ stencil. 

The inequality (1.13) leads to an a'priori pointwise estimate on the 
approximate solution w(x,t). Such pointwise estimates play an essential role 
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with regard to question of the convergence of entropy stable schemes. In 
particular, DiPerna [2, Section 7] has recently shown that in certain cases, 
such (two-sided) estimates are sufficient in order to guarantee the 
convergence of such schemes. 


2. GENERALIZED ENTROPY FUNCTIONS OF THE EULER EQUATIONS 
We consider the Euler equations for polytropic gas 


3 

3t| 


’P 

m 

E 


+ l hr 

k=l k 


m. 


(k) 


q k m + pe 
q k (E + P) 


0 . 


( 2 . 1 ) 


It is well-known, e.g., [1], that for all smooth solutions of (2.1) the 

specific entropy^) 

S(x,t) * lntpp^), 


remains constant along streamlines, i.e., 


DS 3£ v 3S_ 

Dt = 3t + , , q k 3x. 

k=l k 


0 . 


(2.2a) 


Let h(S) be an arbitrary smooth function of S. Multiplying (2.2a) by 
ph'(S) — prime denoting S-differentiation, we find 


(3) 


After normalization, taking the specific heat constant to be c 


1 . 
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P 


3h(S) 

3t 




3h(S) 

3x, 

k 


0 . 


Adding this to the continuity equation which is premultiplied by h(S), 


- 3 3 m, 

H MS) ♦ l £ 


0 , 


( 2 . 2 b) 


we obtain after changing sign, a conservative entropy equation like (1.3) 
which reads [5] 

I 7 [-ph(S)] + l [-t^MS)] - 0. (2.3) 

u k=l k 

In order to comply with the further requirement of being a generalized entropy 

function, U(u) « -ph(S) has to be a convex function of the conservative 

T 

variables u = (p ,»,E) . A straightforward computation carried out by Harten 

[5, Section 2] in the two-dimensional case shows that the Hessian U is 

uu 

positive definite if and only if 


p[h'(S) - Y*h~(s)] > 0. 

Excluding negative densities we may summarize that there exists a family of 
(generalized) entropy pairs (U/£) associated with Euler equations (2*1), 

U(u) = -ph(S), F^^(u) r= -m^h(S) k *= 1,2,3, (2.4a) 

generated by the smooth increasing functions h(S) which satisfy 
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h'(S) - yh"(S) > 0. 


(2.4b) 


3. A MINIMUM ENTROPY PRINCIPLE 

T 

Let u = (p,a,E) be an entropy solution of the gas dynamics equations 

(2.1). Such a solution is characterized by the entropy inequality (1.7) 


9U(u) 2 3F (k) (u) 

9t k=l = 


(3.1) 


which holds for all entropy pairs (U,F) connected with the equations* 

To derive a minimum principle, we shall make use of an argument due to Lax 
[9, Section 3]* We begin with 


Lemma 3*1: Let u be an entropy solution of the gas dynamics equations 
(2*1), Then for all nonpositive smooth increasing functions h(S) satisfying 
(3* 2b), we have 


/ p(x,t)*h(S(x,t))dx > / p(x,0)*h(S(x,0))dx. (3*3) 

|x|<R l*IP +t “W 

Here q max denotes the maximal speed |q| in the domain * 

Proof : As in [10, Theorem 4*1] we integrate the entropy inequality (3.2a) 

over the truncated cone C = {|x| < R + (t - T ^ # 9 max l® < T ^ t} ; if we let 
(n^,n) denote the unit outward normal, then by Greenes theorem 
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/ ph(S)- 
3 C 


n 0 + 


3 

I 

k-1 


q k n k 


3x > 0. 


(3.4) 


The integrals over the top and bottom surfaces give us the difference between 
the left and right-hand sides in (3.3) and by (3.4) this difference is bounded 
from below by 

ph(S). 

mantle 

The result follows upon showing that the last quantity is nonnegative . 
Indeed, since by assumption -ph(S) > 0, this is the same thing as 

3 

"o + Z VY i 0i 

k~l 

on the mantle we have 


n r 


d 

+ I 

k-1 


Vk 


3x. 


and hence 


(n Q ,n) 


A + 


^max 



n o + 


3 

l q k n k 

k=l 


1 




2 

Snax 


L + ? 

|> 

1 

/ 3 
■la - y 

I’ll \ 

l q max ^ , j 


• v max k-i 

in! J 

\ k-1 x | ) 

f 

/l + q 2 

1 4 1 J 


l max 


> 0 


as asserted. 


The discussion in Lemma 3.1 was restricted to smooth function h(S); by 
passing to the limit, its conclusion (3.3) follows for any nonpositive 
nondecreasing function h(S) satisfying (3.2b), whether smooth or not . 
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To derive the minimum entropy principle, we now make a special choice of 
such function, h(S), given by 


h(S) = Min[ S - S Q ,0] , S Q 


Ess inf S(x,0). 
|x|<R+t*q 


max 


(3.5) 


The nonpositive function h(S) is a nondecreasing concave, one, hence 
admissible by (3.2b), and consequently (3.3) applies 


/ p(x,t)«Min[s(x,t) - S 0 ,0]dx > 
|x|<R 


(3.6) 


/ p(x,0)*Min[ S(x,0) - Sg,0]dx. 

lx I <R+t»q 

Now, by the choice of Sq, the integral on the right of (3.6) vanishes since 
Min[ S( x,0)-Sq , 0] does. The inequality (3.6) then tells us that the integral 
on the left is also nonnegative. But since the integrand on the left is by 
definition nonpositive, this can be the case provided this integrand vanishes 
almost everywhere; that is , we have for almost all x, [x| £ R 


and (1.9) follows. 


S(x,t) > S Q 


Ess inf 
|x*| <R+t»q 


max 


S(x, t=0) 


The minimum entropy principle was deduced from the entropy inequality 
(3.2), which in turn was postulated based on the formal regularization 
introduced in (1.4). In general, other regularizations equally apply; in 



- 11 - 


particular, Euler equations are usually sought as the vanishing viscosity 
limit of the Navier-Stokes equations (here we take for simplicity the one- 
dimensional case)^^ 



+ h 

m 

m 

qm + p 

_E_ 


q(E + p)_ 


r o i 



q 


il 

Bx * 

ia 

Bx 


p I 0. 


(3.7) 


Do the (generalized) entropy inequalities (3.2) remain valid on the basis of 

such limit? To answer this question we first note that if U(u) is any 

entropy function, then thanks to its convexity the mapping u + v = is 

one-to-one, and hence one can make the change of variables u « u(v). Harten 

[5] has shown that such change of variables by each member of the family of 

entropy functions (2.4) puts the viscosity terms on the right of (3.7) into a 

negative semidefinite form. This makes apparent the dissipative effect of 

2 

these viscosity terms. Indeed, if T = c • E _ V2*|<|| denotes the absolute 
temperature, then direct manipulation of (3.7) yields, e.g., [1, Section 63], 
[12, Section 6.10] , 


at" [p h ( s )] + [mh( S) ] 


p-h(S) 



(3.8) 


from which we recover the entropy inequality (3.2a) for all smooth increasing 
functions h(S). We note that the convexity condition was not assumed in this 


^^With p combining the two viscosity coefficients in the general Navier- 
Stokes equations. 
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case* The merit of using the convexity condition, however, is that it enables 
us to deal with more general artificial viscosity terras, other than those 
appearing in the Navier-Stokes equations* Such artificial viscosity terras are 
frequently encountered in finite-difference approximations to the Euler 
equations; a specific example of this kind is studied in the next section. 

Finally we would like to remark on the previously mentioned Navier-Stokes 
equations. Our discussion above took into account only the viscosity 
contribution, neglecting heat conduction. Hughes, et al., [7] have shown that 
when the heat flux is also added, compare (3.7), 



"P" 


m 


■ 0 " 


0 


3 

dt 

m 


qm + p 

3 

" V 37 

is. 

3x 

+ K h 

0 

(3.9) 




q(E + p) 

! 

3 q 

L q £ J 


3T 

8x 



with k denoting the heat conductivity constant, then only the ^physical" 
entropy, U(u) *= -pS survives as the one which puts the additional heat flux 
into a symmetric negative-definite form. We would like to note in this 
connection the difference limit behavior of the Navier-Stokes flows depending 
on the viscosity and heat conductivity; Gilbarg [4] has shown that as k + 0 
keeping p fixed, we are led to a continuous thermally nonconducting shock 
layer, whereas for p + 0 with k fixed the convergence is to a (generally) 
discontinuous nonviscous shock layer. Consequently, the viscosity rather than 
the heat flux should play the major rule in an appropriate regularization 
model for the Euler equations. 
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4. DISCRETE APPROXIMATIONS OF THE EULER EQUATIONS 

In this section we consider approximate solutions of the Euler 
equations, w(x^,t), whose total entropy decreases in time, compare (1.10) 

£ v U(w(x v ,t ■+ At))Ax v < ]> v U(w(x v ,t))Ax v . (4.1) 

Estimate (4.1) holds for all entropy functions U ® -ph(s) in (2.4). By 
passing to the limit, this applies to our previous choice of the function 
h(s) in (3.5) 

h(s) = Min[S - S Q ,0] , (4.2a) 

this time with a constant Sq which is taken to be 

Sq = Min S(w(x v ,t)). (4.2b) 

v 

By our choice of Sq, we have U(w(x v ,t)) = 0. The inequality (4.1) tells us 
that the left-hand side is therefore, nonnegative; consequently 

S(x,t + At) - S Q 2 h(S(x,t + At)) 2 0 


and (1.13) follows. 

Approximate solutions which fulfill the required estimate (4.1) can be 
obtained by entropy stable schemes satisfying the cell entropy inequality 
( 1 . 11 ) 


U(w(x v ,t 


+ At)) < U0 


r(x v ,t)) 


d 

I 

k=l 


_L [ F (k) 


Ax 


V+V2 




(4.3) 
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Examples of such entropy stable schemes include the Godunov-type and Lax- 
Friedrichs schemes, e.g., [6]. A more precise minimum principle follows in 
these cases, taking into account the support of the schemes' stencil. In 
particular, the (one-dimensional) Godunov scheme results from averaging of two 
neighboring Riemann problems [6], each of which satisfies (1.9). Consequently 
we have the 


Minimum Principle (of the Godunov scheme): Let w(x^,t) the Godunov 

approximate solution to the Euler equations (2.1). Assume that the 

appropriate CFL condition is met. Then the following estimate holds 

S(w(x ,t + At)) > Min S(w(x.,t)). (4.4) 

v-Kj<v+l J 

Since the Lax-Friedrichs scheme coincides with a staggered Godunovas solver, 
the same conclusion, (4,4), holds. Another way to see this is outlined below; 
it makes no reference to Riemann^s solution and can be generalized to the 
multidimensional problem. 

To this end, we approximate the (for simplicity — one-dimensional) Euler 
equations with the Lax-Friedrichs scheme 


w(x ,t + At) 
v 


M w(x v+i ,t) +w( vi’ t) ] 

■7 " fO^v-i’^) > 



(4.5) 


We remark that the Lax-Friedrichs scheme can be derived from center 
differencing of the regularization model (1.4) • Lax has shown [9, Theorem 
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1.2] that if X = ~ is sufficiently small, then solutions of this difference 
scheme satisfy the following cell entropy inequality 


U(w(x v ,t + At)) < 


U(w(x v+1 ,t)) + U(w(x v _ 1 ,t)) 
2 


- i [F(w(x v+1 ,t)) - Wx v _ lt t))] 


(4.5)- 


for all entropy pairs (U,F) « (-ph(S) ,-mh(S)) in (2.4). by passing to the 
limit, this applies to our previous choice of the function h(S) in (3.5) 


h(S) = Min[ S - S Q ,0], 

this time, with a contant Sq which is taken to be 

S 0 = Min[S(x v+1 ,t), S(x v _ 1>t )]. 


(4.7a) 


(4.7b) 


The inequality (4,6) now reads 


p(x v ,t + At)*h(S(x v ,t + At)) > 


1 + Xq(x , ,t) 

2 P( x v _ 1 » t > ,h ( s ( x v _ 1 » t >) 


1 ~ *q( x )+1 »t) 

+ 5 p(x v+1 ,t)-h(S(x i+1 ,t)) 


V + l ! 


(4.8) 


By our choice of the funtion h(S) in (4.7), we have h(S(x^ ± ^,t)) e 0. The 
inequality (4.8) tells us that the left-hand side is therefore nonnegative; 
consequently 
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0 < h(S(x v ,t + At)) < S(x,t + At) - S Q 
and the following minimum principle follows 

S(w(x v ,t + At)) > Min S(w(x v±1 ,t)) . 
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